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Abstract 

In this paper, we study the asymptotic behavior, as the time t goes to zero, of the 
trace of the semigroup of a killed relativistic a-stable process in bounded C 1 ' 1 open 
sets and bounded Lipschitz open sets. More precisely, we establish the asymptotic 
expansion in terms of t of the trace with an error bound of order f 2 / a t- d / a for C 1 ' 1 open 
sets and of order i L l a i~ d l a fo r Lipschitz open sets. Compared with the corresponding 
expansions for stable processes, there are more terms between the orders t~ d l a and 
£(2-d)/a £ Qr p en se t S) aric j ; when a £ (0, 1], between the orders t~ d l a and t^^ d ^ a 
for Lipschitz open sets. 

1 Introduction and statement of the main results 

For any m > and a G (0,2), a relativistic a-stable process X m on R d with mass m is a 
Levy process with characteristic function given by 

E [exp(i£ ■ (X? - X m ))] = exp(-t((|£| 2 + m 2 / a ) a / 2 - m)), £ G R d . (1.1) 

The limiting case X°, corresponding to m = 0, is a (rotationally) symmetric a-stable process 
on W L which we will simply denote as X. The infinitesimal generator of X m is m — (m 2 / Q — 
A)"/ 2 . Note that when m = 1, this infinitesimal generator reduces to 1 — (1 — A) Q / 2 . 
Thus the 1-resolvent kernel of the relativistic a-stable process X 1 on M. d is just the Bessel 
potential kernel. When a = 1, the infinitesimal generator reduces to the so-called free 
relativistic Hamiltonian m — a/— A + m 2 . The operator m — \J — A + m? is very important 
in mathematical physics due to its application to relativistic quantum mechanics. 

In this paper, we will be interested in the asymptotic behavior of the trace of the semi- 
group associated with killed relativistic a-stable processes in open sets of M d . The process 
X m has a transition density p m (t, x, y) = p m (t, y — x) given by the inverse Fourier transform 

P m (t,x) = (2n)- d [ e ~i^ e -t(\e+my^ + mt d ^ 
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For any open set D in M d , the killed relativistic a-stable process X™' 13 is denned by 



m,D 



X? if t<T#, 

d if t > r™, 



where rg 1 = inf{t > : X t m ^ .D} is the first exit time of X m from D. The process X t m ' D is 
a strong Markov process with a transition density p^(t,x,y) given by 



pZ(t,x,y)= P m (t,x,y)-r%(t,x,y), 



with 



rS(f , x, y) = E x [t > r£; p m (t - r£, X™,y) . 
We denote by {P™' D : t > 0) the semigroup of Xf on L 2 (D): for any / G L 2 (D), 



ym,D 



f(x) := E, f{X? D ) 



f{y)p%{t,x,y)dy. 



D 



Whenever D is of finite volume, P™' D is a Hilbert-Schmidt operator mapping L 2 (D) into 
L°°(D) for every £ > 0. By general operator theory, there exist an orthonormal basis of 
eigenf unctions {4>^}^ =1 for L 2 (D) and corresponding eigenvalues {An }J£=i of the generator 



of the semigroup P^' D satisfying 



< AS m) < } < A< m) < 



with An'™' 1 — > oo. By definition, we have 



We also have 



An will be simply denoted by A n . 

In the remainder of this paper, we assume d > 2. We are interested in finding the 
asymptotic behavior, as t — > 0, of the trace defined by 

« oo « oo 

Zn(t)= / ^(t,*,*)cix = Ve^' / (^ m) ) 2 (^ = E e ~ A ^ 

jD n=l jD n=l 

It is shown in [2] that for any open set D of finite volume, it holds that 



Pd 



{t^y) = Y,e~ X ^^\x)<P { r\y)- 



x £D,t>0. 



n=l 



limt d/a Z° D = C 1 \D\, Ci 

t->-o 



(2Tt) d a 



;i-2) 



This is closely related to the growth of the eigenvalues of P 4 °' D : if N°(X) is the number 
of eigenvalues Xj such that Xj < X, then it follows from the classical Karamata Tauberian 
theorem (see for example pi)]) that 

This is the analogue for killed stable processes of the celebrated Weyl's asymptotic formula 
for the eigenvalues of the Dirichlet Laplacian. We will see later in this paper that exactly the 
same formula is true for relativistic stable processes. That is, the first term in the expansion 
of Zp(i) is the same as that of Zp(t) and (jl.3p is also true for relativistic stable processes. 

Our main goal in this paper is to get the asymptotic expansion of Z™(t) as t — > 
under some additional assumptions on the smoothness of the boundary of D. Our work is 
inspired by the paper [7] for Brownian motion and the papers [21 [3] for stable processes. 
The first theorem is an asymptotic expansion of Z™(t) with error bound of order t 2 / a t~ d / a 
in C 1,1 open sets. To state the result precisely, we need some definitions. Recall that an 
open set D in M, d is said to be a (uniform) C 1 ' 1 open set if there are (localization radius) 
R > and A such that for every z G dD, there exist a C 1 ' 1 function <fi = <p z : W d — > 
R satisfying 0(0, •• • ,0) = 0, V0(O) = (0, ...,0), |V0(x) - V(j)(y)\ < A \x - z\ and an 
orthonormal coordinate system CS Z : y = (yi, ■ ■ ■ ,yd-i,yd) '■= {ViVd) with origin at z such 
that B(z, R)f)D = {y = (y, y d ) G 5(0, R) in CS Z : y d > <j)(y)}. For x G M d , let 6 D (x) denote 
the Euclidean distance between x and D c and 5qd{x) the Euclidean distance between x and 
dD. It is well known that a C 1,1 open set D satisfies both the uniform interior ball condition 
and the uniform exterior ball condition: there exists r$ < R such that for every x G D with 
o~8d(x) < r o an d y G M. d \D with Sgpty) < r , there are z x ,z y G dD so that \x — z x \ = Sqd{x), 
\y — z y \ = SoD^y) and that 5(xo,r ) C D and B(y ,r ) cM d \D, where Xq = z x + r (x — 
z x )/\x — z x \ and yo — z y + r (y — z y )/\y — z y \. In fact, D is a C 1 ' 1 open set if and only if D 
satisfies the uniform interior ball condition and the uniform exterior ball condition (see [H 
Lemma 2.2]). In this paper we call the pair (ro,Ao) the characteristics of the C 1 ' 1 open set 
D. For any open set D in M. d , we use \D\ to denote the d- dimensional Lebesgue measure of 
D and ^^(dD) to denote the (d — l)-dimensional Hausdorff measure of dD. When D is 
a C 1 ' 1 open set, H a!_1 ((9-D) is equal to the surface measure \dD\ of dD. We will use H to 
denote the half space {x = (x±, X2, ■ ■ ■ , Xd) '■ x\ > 0}. 

The following is the the first main result of this paper. 

Theorem 1.1 Suppose that D is a bounded C 1,1 open set in M. d . Let k be the largest integer 
such that k < -. Then the trace Zp(t) admits the following expansion 

7^(f\ r\n\i-~ r\f>r>\+— . u d T(d/a)\D\ £ A m n n( t 2 ' a 
Z D (t) = C l \D\t o-C 2 \dD\t« + a t "2^1^ + ^)' 

^ ' 71=1 
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where C\ is given in (jl.2p and 



Co 



(r,6), (r,0))dr. 



The second main result of the paper is an asymptotic expansion of Z^it) with error 
bound of order t x l a t~ d l a i n Lipschitz open sets. Before we state the second main result, we 
recall the definition of Lipschitz open sets. An open set D in ~R d is called a Lipschitz open 
set if there exist constants Rq (localization radius) and A > (Lipschitz constant) such that 



for every z e dD there exist a Lipschitz function F : 



Dd— 1 



with Lipschitz constant A 



and an orthornormal coordinate system y = (yi, ■ ■ ■ ,y d ) such that D fl B(z, Rq) = {y : y d > 
F{yii " " • > Vd~i)} H B(z, Rq). Here is the second main result. 

Theorem 1.2 Suppose that D is a bounded Lipschitz open set in ~R d . Let j be the largest 
integer such that j < -. Then the trace Z^{t) admits the following expansion 



t d/a Z%(t) = d\D\ - C 2 H d ~\dD)t 1/a + 



co d T(d/a)\D\ 
(2ir) d a 



n=l 



where C\ and C 2 are the same as in Theorem 

The asymptotic behaviors of the trace Z^{t) of the killed Brownian motion (i.e., killed 
symmetric a-stable process with a = 2) in bounded domains D of M. d have been extensively 
studied by many authors. It is shown in [5] that, when D is a bounded C ' domain, 

c\D\t l - d ' 2 



Z D {t) - (4nt)- d / 2 
The following asymptotic result 

Z D (t) = (Ant)- d / 2 



D 



\dD\ 



< 



R 2 



t > 0. 



D\ 



\dD\ + o(t 



1/2- 



t 0, 



was proved in [6] when D is a bounded C 1 domain. (jl.4p was subsequently extended to 
Lipschitz domains in [7]. 

The asymptotic behaviors of the trace Z° D (t) of killed symmetric a-stable processes, 
< a < 2, in open sets of M. d have been studied in [21 E]. It was shown in [2] that, for any 
bounded C 1 * 1 open set D, 



Z° D (t) 



d\D\ C 2 \dD\t^ a 



t d/c 



+ 



t d/c 



< 



\D\t 2 / c 

r 2 t d/a 



where C\ and C 2 are the same as in Theorem 1 1 . 1 1 and c is a positive constant depending on d 
and a only. It was shown in [3] that, when D is a bounded Lipschitz domain, Z° D (t) satisfies 



t d/a Z° D (t) =C X \D\ 



C 2 H d -\dD)t l/a + o(t 1/a ) 



4 



Remark 1.3 Note that the first term in the expansion of Z^(t) is exactly the same as in 
the case of Z° D {t). However the rest of the terms are quite different. We note here that the 
coefficient of the term of order t 1 / a t~ d / a is the same in the stable process case, but in the 
case of relativistic stable processes for C 1,1, open sets, there are k intermediate terms of the 
form t k t~ d l a , where k is a positive integer less than 2/ a. Since < a < 2, there is at least 
one more term involved in the asymptotic expansion of Zp(t) than that of Z° D (t) up to order 
oft 2 l a t- d l a . ForLipschitz open sets, when a < 1 there are j intermediate terms of the form 
tH~ d l a , where j is an integer that is less than or equal to 1/a. 

Remark 1.4 In an asymptotic expansion for the trace of relativisitic a-stable processes 
in bounded C 1,1 open sets was established. Compared with Theorem \l.l[ the expansion of J^j 
does not contain the intermediate terms. 

The rest of the paper is organized as follows. In Section O we recall some basic facts 
about relativistic stable processes and present several preliminary results which will be used 
in Sections [3] and HI Theorem 11.11 is proved in Section [3j while Theorem 11.21 is proved in 
Section HI 

Throughout this paper, we will use c to denote a positive constant depending (unless 
otherwise explicitly stated) only on d and a but whose value may change from line to line, 
even within a single line. In this paper, the big O notation f(t) = 0(g(t)) always means 
that there exist constants C and to > such that f(t) < Cg(t) for all < t < to- 

2 Preliminaries 

In this section, we recall some basic facts about relativistic a-stable processes. From (II. ip . 
one can easily see that X m has the following approximate scaling property: 

{m- 1/a (X^ t - Xq 1 ), t > 0} has the same law as {X™ - X m , t > 0}. 

In terms of transition densities, this approximate scaling property can be written as 

p m (t, x, y) = m d/a p\mt, m 1/a x, m l/a y). (2.1) 

It is well known that the transition density p"^(t, x, y) of X m,D is continuous on (0, oo) x 
D x D. Since both p m (t, x, y) and p^{t, x, y) are continuous on (0, oo) x D x D, r^{t, x, y) = 
p m (t,x,y) — p^it.x.y) is also continuous there. pp(t,x,y) and r^{t,x,y) also enjoy the 
following approximate scaling property: 

p l ml/aD (t, x, y) = m~ d ' a p™(t/m, x/m 1 ^, y/m 1 ^), (2.2) 

r l m i /aD (t, x, y) = m-^r^t/m, x/m 1 ^, ylm l l a ). (2.3) 
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The Levy measure of the relativistic ct-stable process X m has a density 

POO 

J"(s) = n\x\) := a / ^u)-"l\-^^e-^u-^d^ 
21 (1 - a/2) J 

which is continuous and radially decreasing on M d \ {0} (see [T3J Lemma 2]). Put J m (x, y) : = 
j m (\x -y\). Let A(d, -a) := a2 a - 1 Tz~ d / 2 Y{^)V{\ - f)" 1 . Using change of variables twice, 
first with u = \x\ 2 v then with v — 1/s, we get 

J m (x, y) = A(d, -a)\x - y\- d -^{m l ' a \x - y\), (2.4) 

where 

V(r) := ( — ^ 1 /°° S ( rf+ ")/ 2 - 1 e- s / 4 -'- 2 /^ S , (2.5) 



which satisfies ^(0) = 1 and 

c -l e -r r (d+a-l)/2 < ^ < Cie - v (d +a -l)/2 Qn ^ ^ 

for some ci > 1 (see [9j pp. 276-277] for details). We denote the Levy density of X by 

J(x, y) := J°(x, y) = A(d, -a) \x - y\~ d ~ a . 
Note that from (|2~4"1) and ( 1231) we see that for any x <ER d \ {0} 

r(N)<j°(ix|). 

It follows from [HI Theorem 4.1.] that, for any positive constants M and T there exists a 
constant c > 1 such that for all m G (0, M], t G (0, T], and x, y G M d we have 

c -l ^-d/a A t jm^ y ^ < p m^ X; y ) < c ( r «*/« A £ /"(^ y )) . ( 2 .6) 

We will need a simple lemma from [TT] about the relationship between r^(t,x,y) and 
r%(t,x,y). The lemma is true in much more general situations but we just need it when 
one of the processes is a symmetric a-stable process and the other is a relativistic a-stable 
process. 

Lemma 2.1 Suppose that X and Y are two Levy processes with Levy densities J x and J Y , 
respectively. Suppose that a — J x — J Y is nonnegative on W 1 with f Rd a(x)dx = £ < oo and 
D is an open set. Then for any x G D and t > 0, 

Pd&x,-) < e a p^{t,x,-) a.s. 

If, in addition, p x (t, ■) and p Y (t, ■) are continuous, then we have for x,y G D, 

rl{t,x,y)<e m r x {t,x,y). 
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The next proposition is the (generalized) Ikeda-Watanabe formula for the relativistic 
stable process, which describes the joint distribution of and X r T 



m 
T D 



Proposition 2.2 (Proposition 2.7 [12J) Assume that D is an open subset of M. d and A 
is a Borel set such that A C D c \ dD. If < t\ < t2 < oo, then 



¥ x eA, tx < t% < t 2 ) = J J* p"%(s,x,y)ds J J m (y,z)dzdy, x G D. 

Now we state a simple lemma about the upper bound of r^(t, x, y), which is an analogue 
of [21 Lemma 2.1] for stable processes. 

Lemma 2.3 Let M,T be positive constants. Then there exists a constant c = c(d,a, M,T) 
such that for all m G (0, M] and t G (0, T] we have 

r-(t xvXc (r^ A *^™ 1/a M*)) 
r D {t,x,y)<c\t A 

Proof. Since if) is eventually decreasing and ^(0) = 1 > 0, there exists a constant c\ > 
such that ip(x) < c\ip(y) for all < y < x. Now from the definition of r^(t, x, y) and (I2.6P 
we have 

r%(t,x,y) = r%(t,y,x) 



< E„ 



< E„ 



t>T^ ]P m (t-T^,X^,x) 

til){m 1/a \x - X™ m \) 



c t 



-d/a 



A 



\X 



ym \d+a 



< CCi [t 



d/a a f#ng^ 



S D (x) 



d+a 



□ 



We will need two results from [2]. The first result is about the difference pp(t,x,y) — 
p^(t, x, y) when D C F. The proof in [2J, given for stable processes, mainly uses the strong 
Markov property and it works for all strong Markov processes with transition densities. 

Proposition 2.4 (Proposition 2.3 |2J) Let D and F be open sets in M. d such that D C F. 
Then for any x, y G M. d we have 



p™{t, x, y) - pS(t, x, y) = E, r™ < t, X™ G F \ D : p™(t - X%,y) 



7 



Now we introduce some notation. Recall that if D is a C ' l open set with characteristics 
(ro, Ao), then for every x G D with Sqd{x) < r and y G M, d \ D with Soniy) < r , there 
are z x ,z y G 3D so that \x — z x \ = Sqd(x), \y — z y \ = Sgr>(y) and that B(xo,r ) C D and 
B(y ,r ) C lR d \Z), where x = z x +r (x-z x )/\x-z x \ and y = ^+ r o(2/-%)/|y-^|- Let H(x) 
be the half-space containing B(xo,r ) such that dH(x) contains z x and is perpendicular to 
the segment z x z y . The next proposition says that, in case of the symmetric a-stable process, 
for small t, the quantity r^(t, x, x) can be replaced by r° H ^{t, x, x), which was a very crucial 
step in proving the main result in [2]. 

Proposition 2.5 (Proposition 3.1 of |2J) Let D C M. d be a C 1 ' 1 open set with charac- 
teristics (r ,A ). Then, for any x with 5qd(x) < r /2 and t > with t l / a < r /2, we 
have 

cf l/a ( / t l/a 



m,x,x)-^ x) (t,x,x)\<^ [{s^j) Al 

We will need some facts about the "stability" of the surface area of the boundary of C 1,1 
open sets. The following lemma is [HI Lemma 5]. 

Lemma 2.6 Let D be a bounded C 1,1 open set in M. d with characteristic (ro, Ao) and define 
for < q < r Q , 

D q = {x G D : 8 D (x) > q}. 

Then 

)d—l / \ d— 1 

\dD\<\dD q \< (-^-j 0<g<r . 

The following result is j2j Corollary 2.14]. 

Lemma 2.7 Let D be a bounded C 1,1 open set in M d with characteristic (r ,A ). For any 
< q < r /2, we have 



1. 2~ d+l \dD\ < \dD q \ < 2 d - 1 \dD\, 

2. \dD\ < 



2 d \D\ 



3. \\dD q \-\dD\\<*^<*^\. 
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3 Proof of Theorem 11.1 



We first prove that lim^o t a Z^{t) exists and identify the limit. 

Lemma 3.1 The limit lim^o t : " Z^(t) exists and is equal to Ci\D\, where C\ is the constant 
in Theorem ! 

Proof. By definition, 

t d/a Z™{t) = t d/a [ p%{t,x,x)dx 
= t d/a ( [ p m (t,x,x)dx- [ r%(t,x,x)dx J . (3.1) 



D J D 



For the first integral on the right hand side of (13. ip . note that, by the approximate scaling 
property (12.21) and the dominated convergence theorem, we have, as t — > 0, 

t d/a Q p m (t,x,x)dx\ = J p tm {l,x,x)dx =\D\p tm (l,0) 

It remains to show that lim^o t d ^ a j D r 1 £(t,x,x)dx = 0. By Lemma [2.31 we have that 
t d ' a r™{t, x, y) < c, (t, x, y) 6(0,1] x fix D, 
for some c > 0. Hence we have by the monotone convergence theorem, 

t d/a / r%(t,x,x) ^ as t -> 0. 

JD\D tl/2a 

For x G D t i/2 a we have by Lemma 12731 again for t G (0, 1], 

r™(t,x,x) < ct^ + ^ , x G D t i/ 2a . 
Hence lim t ^ t d l a J D ^ r^(t, x, x)dx = 0. □ 

It follows from Lemma 13.11 that if N m (X) denotes the number of eigenvalues X^ such 
that X" 1 < A, then it follows from the classical Karamata Tauberian theorem (see for example 
[ID]) that 

N m (X) ~ ^fl, ^ X d/a , as A oo. 
v ' T{d/a + l) 
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This is the analogue for killed relativistic stable processes of the celebrated Weyl's asymptotic 
formula for the eigenvalues of the Dirichlet Laplacian and it is already proved in [4] (see [U 
(1.10)]). This result has been known at least since 2009, see jU Remark 1.2]. 

Now we focus on identifying the next terms in Z^(t). For this, we need to find the order 
of t in Z%(t) - CiT*. Note that by Lemma EEJ 

Z£{t) - dr dla = [ p%{t, x, x) - p°{t, x, x)dx 
Jd 

= / (p m (t, x,x) —p°(t, x,x)) dx - / rp(t,x,x)dx. 
Jd Jd 

The next lemma gives the orders of t in p m (t, x, x) — p°(t, x, x) up to t"t"«. 
Lemma 3.2 Let k be the largest integer such that k < —. Then we have 

P m (t, X, X) - P % X, X) = t -d/ a ^r( d M y H^ t n + 0{ j2/a t -d/ a y 

[2it) d a> ' nl 

v ' n=l 

Proof. By the scaling property (12.11) we have 

p m {t, x, x) - p°(t, x, x) = p m {t, 0) - p°{t, 0) 

= r d/a (p tm (i,o)-/(i,o)) 



Note that for any x > we have (1 + x) a ^ 2 < 1 + fx. Thus 

(i £ i> +(mt rr=ifi-(i + ¥)^<i«i-(i + ^ 



Consequently 



< e~ |J;l -e 



ier _ P -(iei 2 +M) 2/< T /2 



< e 



-l?l c 



a (mtf /a ' 



2 iei 2 -° y 

where we used 1 — e _2; < x for all x > in the last inequality above. Therefore 
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< 



< 



< 



e -(\ti\ 2 + (™t) 2/a ) a/2 +rnt _ e -\t\ a e mt + g-|€|° e mt _ e ~\t\ a 
e -(l€| 2 + M) 2/a ) Q/2 +mt _ e -\t\ a e ™t 



# + / | e -ier e mt -e-i«i Q i 
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lei 



2-Q 



n=l 



(mt) ? 



n=fc+l 



Since + j > 2/a for any j > 1 , we have 

-(\£\ 2 +(mt) 2 / a } a/2 +mt 



(mi)' 



rv. 



0(t 2/a ). Therefore 



- e-w) di = o(tv*) + UdT{d/a) y ^ 

/ a ^— ' n! 

' n=l 



and 



p m (t, x, x) - p°(t, x, x) = t - d '" WdTWa) 



(27r) d a 



J2^ t n + ^ t -d/ay 



n=l 



□ 



Now we try to find the orders of t in the expansion of J D r^(t,x,x)dx up to the order 
of £«t~a. For this, we need to assume some regularity condition on the boundary of D. 
Hence in the remainder of this section we assume that D is a bounded G ,x open set with 
characteristic (r ,A ). We also assume that t l l a < y. 

We first deal with the contribution in T) ro i%. 

Lemma 3.3 There exists c = c(d, a) > such that 

[ r£(t,x,x)dx<ce 2 -M^. 

Proof. It follows from Lemma I2TT1 that r™(t,x,y) < e 2mt r® D {t,x,y). By [2, (3.2)] we know 
that 

c\D\t 2 / a 



r° D (t,x,y)dx< . 



r /2 



The desired assertion follows immediately 



□ 



Lemma 3.4 There exists c = c(d, a) > such that 



rp(t,x,x) r H(x)(tj x i x ) — ce2mt ^dJ 



5 D (x) 



d+%-1 



A 1 
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and 



/ (r£(t, x, x) - r™ {x) (t, x, x)) dx < 

JD\D rn/ o 



ce 



2mt 



t 2/a 



'D\D rQ/2 

Proof. If the first assertion of the lemma is right, then it is easy to see that 



L 



t l/a \ 



M]dx< ct 1/a . 



>D\D ro/2 \\$d(x). 

Hence we focus on proving the first assertion. By [2j (3.4)], we know that 



r%(t, x, x) — r H ^ (t, x, x) < c 



t d ' a \ \5 D {x) 



Al . 



Recall that J m (x) < J°(x) for any x G M. d \ {0}. Now it follows from the generalized 
Ikeda-Watanabe formula and Lemma 12.11 that 

r D (t, x, x) — f H ^ (t, x, x) 



E, 



* > r™, X% G H(x) \ D; p™ [x) (t - t%, X™ 



x 




Po(s,x,y)ds / J m (y,z)PH( x )(t - s,z,x)dzdy 

'/>./(> JH(x)\D 

< e 2mt I I p° D (s,x,y)ds I J°(y, z)p° H{x) (t - s, z, x)dzdy 




D JO 



e 2m % 

2mt (JO 



JH(x)\D 

t > r% , X T o G H(x) \ D; p° H(x) (t - r° D , X T% , 



X 



r° D (t, x, x) — r° H ^(t, x, x)) 
< r p 2mtt — ( (J——\ d +^- 1 A ]\ 



Lemma 3.5 There exists c = c(d, a) > such that 



□ 



r% {x) (t,x,x)dx - t 1/a r d/a / 24 \dD\f^{l,u)du < ct 2 ' a t- d ' a . 



'D\D rQ/2 

Proof. Using the scaling relation (12. 3 p we get 



J D\D,, o/2 
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r /2 



ro/2 



\dD u \f™(t,u)du 

\dD u \r d / a f^(i,u/t^ a )du 







rr /2t 1 / a 

t 1/a t~ d/a / \dD utl/a \f^(l,u)du. 
Jo 



It follows from CorollaryOthat \ \8D q \ - \dD\ \ < 



2 d dq\dD\ ^ 2 2d dq\D 



D\D rQ/2 



< t 1/a r d/c 



r H(x) (A x i x ) ~ t 1 a t d / 



'•q 

" ' 2(1 " 19^1/^(1, u)du 



for any g < Tq/2. Hence 



\\dD utiJa \-\dD\\f™{l,u)du 



< Cl t 2/a t- d/a / ufjF(l,u)du 

Jo 

< c 2 t 2/a t- d/a . 



Lemma 3.6 There exists c = c(d, a) > such that 

t l/a t -d/a / \dD\f^(l,u)du _ t l/* t -d/* / 2t " \dD\f t ^ l (l,u)du < ct 2 '^^ 

Jo Jo 
Proof. It follows from Lemma [2. II that 

rn 

jX/a^—- d/a 



\dD\f™(l,u)du-t 1/a r d/a I ^ \dD\f^(l,u)du 
'o Jo 

/EX) 

2t 1 /" 

/OO 

2t l/a 

/OO 
f H {l,u)du. 

2t 1 '<* 

For q > r /(2t 1/a ) we have g) < cq d a < cq 2 . Hence 



/// 



f H {l,u)du<c -|<c 



2t 1 / a 



r ? 

2t 1 / a 



r 



and the result now follows. 



□ 



□ 
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/OO /'DC 
f t ^(l,u)du = J f H {l, 



u)du. 



Proof. This follows immediately from the continuity of m t- > r^(t, x, y) and the dominated 
convergence theorem. □ 



Proof of Theorem 11.11 Combining Lemmas 13.11 I3.2[ 13.31 13. 4[ 13. 5[ 13. 6[ and 13. 7[ we 

□ 



immediately arrive at Theorem 11.11 



4 Proof of Theorem 1.2 



In this section we always assume that D is a bounded Lipschitz open set in M. d . The argument 
of this section is similar to previous section and [3] . We will follow the argument in [3] closely, 
making necessary modifications for relativistic stable processes. Note that even though the 
main theorem in [3] is stated for a Lipschitz domain, it remains true for a bounded Lipschitz 
open set. 

First we need two technical facts which play crucial roles later. The first proposition is 
[31 Proposition 2.9] and we will state it here for reader's convenience. 

Proposition 4.1 (Proposition 2.9. [3]) Suppose that f : (0, oo) — > [0, oo) is continuous 
and satisfies f(r) < c(l Ar" 13 ) for some (3 > 1. Furthermore, suppose that for any < R\ < 
R 2 < oo, f is Lipschitz on [Ri,R 2 ]. Then we have 

limiy f(^^J dx = H d - 1 (dD) J™ f(r)dr. 

Lemma 4.2 Suppose that f : (0, oo) — > [0, oo) is continuous and satisfies f(r) < Ci(l Ar _/3 ) 
for some (3 > 1. Furthermore, suppose that for any < R\ < R 2 < oo, f is Lipschitz on 
[R±, R 2 ]- Let {f v : rj > 0} be continuous functions from (0, oo) to [0, oo) such that, for any 
< L < M < oo, lim/ r? (r) = f(r) uniformly for r G [L, Ml. Suppose that there exists 

c 2 > such that f v (r) < c 2 f(r) for all rj < 1. Then we have 

lim - [ p( M^I] dx = H d -\dD) r f(r)dr. 
v->o+ V Jd \ V J Jo 

Proof. Let i^ v (r) = rj^ 1 \{x G D : Sd{x) < rjr}\. Note (cf. proof of [TJ Proposition 1.1]) that 
ipT]{ r ) < c for all r), r > and that 
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and 



T)- 1 / P 

'd 



7] 



dx 



P(r)d^(r). 



It was shown in [3, Proposition 2.9.] that, for any < Ri < R 2 < oo and r\ > 0, / satisfies 

-Hi 

(4.1) 

, ; -i ,.r> 1 -'- (4_2) 



f(r)diP v (r) < cRx, 

POO 

/ f(r)dip v (r) < cr/ 3 " 1 + cR\~ 

JRi 



lim 



H 2 

R 2 



f(r)di; v (r) = H d - 1 (dD) I f{r)dr. 



Ri 



R-2 



Ri 



Since f v < C2/ for rj < 1 we have the same inequalities as (14.11) and H4.2[) for J' 7 , 77 < 1. 
Hence it is enough to show that 



lim 

^0+ JRx 



Ri 



f v (r)dip, n (r) =U d -\dD) I f(r)dr 



R 2 



Ri 



For any partition R\ = xq < x\ < ■ ■ ■ < x n = R2 of [R\, R2], we have 

n n 



i=l 

n 



i=l 



i=l 



< II P -f lU-ifliAi ^(-^2)- 

Note that for any 77 > the function r — > ip v {r) is nondecreasing and for any 77 > 0, r > 
we have tp v {r) < cr for some constant c. Since — » / uniformly on r G [i?i,i?2], taking 
supremum for all possible partitions gives 



lim 

r/-S-0+ 



H 2 



Rx 



T?->0 + 



H 2 



P(r)d^ v {r)= lim / f(r)dip ri (r) = H d - l (dD)l f{r)dr. 



Hi 



H 2 



Hi 



□ 



Lemma 4.3 For any < L < M < oo ; p m (t,x,y) converges uniformly to p°(t,x,y) as 
m ^ for (t, x, y) G [L, M] x 



Proof. Note that 

\p m (t,x,y) -p°(t,x,y)\ 



(27T) 



fl -tf(i/-<c) /-i((|£p +m 2/o<)<»/2_ m ) _ e ~*|?|' 
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< (2ixy d I e~ i£ - {v ~ x) (V' ((|5|2+m2/Q)a/2_m) _ e -* |?|aN ) d£ 

~ JR d ^ ' 

< (2n)- d [ e -m\ 2 +m^)^- m ) _ e -m^ 

= (2n)- d (p m (t,0)-p°(t,0)). 
Now it follows from Lemma [3.21 that for t G [L, M] and x, y G M. d , 
\p°(t,x,y) -p m (t,x,y)\ 

A JR d Isl ~y n - JR d 

/ JR d Is I n=1 n! JRd 

The last quantity above converges to as m — > 0. □ 
For convenience, we define the following notation. 

/2(f,r):=r£(*,(r,0),(r,0)), r > 0. 
Lemma 4.4 For any < L < M < 00 and m > 0, 

lim/^ n (l,r) = fjj(l, r), uniformly in r G [L,M], 
i/iai is, given e > t/iere exists to > snc/i t/iat /or < t < to we fta^e 

sup |r^(l, (r, 0), (r, 6)) - r° H (l, (r, 0), (r, 0)) | < e. 

r6[L,M] 

Proof. Recall that r° H {t,x,y) = E x [t^ < t,p°(t - r°,X T £,w)] and r™{t,x,y) = E x [t% < 
t,p m {t - T%,X™m,y)}. It is well known that 



p\t,x,y)^r d / a A 



\x - y\ d + a ' 

Since \X® — (r, 0) | > L, we have, together with Lemma [2.11 

p»(l-r«,A-»„,(r,0))<c^, 

1 _ _im 

ptm{1 _ ^xfy, (r,0)) < ce im -^. 
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Now take S\ small so that 



E 



(r,6) 



l-5 1 <r H <l,p (l-T H ,X° T o,(r,0)) 



< e, 



E 



(r,0) 



1 - Si < rfr < l,p tm (l - t^,X^, (r,0)) 



< e. 



(4.3) 



(4.4) 



Now let V m be a Levy process with Levy density a = J — J m and define T m := inf{t > 
: V™ 7^ 0}. Then V m is a compound Poisson process and T m is an exponential random 
variable with parameter m and independent of X (See [13]). Then we have 



E 



(r,0) 
+ E (r,0) 



r*T < 1 - 5x^^(1 - r*T, (r, 0)) 

T tm <^ T t™ <i_ 5 1: p tm (l - T fP, X*^, (r, 0)) 
Since p tm (l - rg 1 , X%, (r, 0)) < c-^, we have 



E 



(r,6) 



r tm < i jT ^i < i _ 5 1?p tm (l - r £», X^, (r, 0)) 



< C 



-(l-e- mt ). 



Similarly we also have 



%,o) [T tm < 1, 4 < 1 - - r° , X° TH , (r, 0))] < c^(l - e 



-mt\ 



(4.5) 



(4.6) 



Take ti > such that (14. 5p and H4.6[) is less than e for all t < t\. Next note that for 
T tm > 1 and rjp < 1, we have r^™ = r£ and X*^ = X° . Hence it follows that 



l%,Q) 

-E 



(r,0) 



< E 



(no) 



T* m > 1, rj? 1 < 1 - 5i,p* m (l - r*T, X^, (r, 0)) 
V m > 1, r° H < 1 - 5 1; p (l - r° H , X°, , (r, 6))] | 
T te > 1, r° < 1 - fc, \p tm (l - r° , X T V (r, 0)) - p°(l - r° , X° , , (r, 0))| 



< sup |p* m (s,a;,?/) -p (s,x,y)|. 

se[5i,l],a;,yeIR d 



(4.7) 



It follows from Lemma 14731 that there exists t 2 > such that sup ae r,y ji x )2/eR <* |p* m (s, x, £/) — 
p°(s, x, y)\ < e for < t < t 2 . Now let t = t± A t 2 . Then for any < t < to we have from 

(ELS), gaD, gSD, gSJ), and gZD 

|r*T(l, (r, 0), (r, 0)) - ^(1, (r, 0), (r, 0))| 
= l%,o)Hr < l,p fm (l-T!r,XX(r,0))] -E (r>5) [r° < l,p°(l - r° , X r °o , (r, 0))]| 
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< |E (rtB) [l > r|T > 1 - 5 1: r^ < l,p tm (l - r^,Xlf m , (r,0))]| + 
|E(r,5)[l >t° h >1-S 1 ,t° h < l,p°(l - r°,X r °o , (r,0))]| 



IE 



(r,0) 



T im < 1, rg 1 < 1 - 5i,p* m (l - r|T, (r, 0)) 
T tm < 1, r° H < 1 - 5 1)P °(1 - r° , X°, , (r, 0))J | 
T tm > 1, r° < 1 - - t%, X% , (r, 0)) 



IE, 



+l%,6) 
-E 



+ 



T im > 1, t£ < 1 - 5 1; p u (l - r° H , X T ° , (r, 0)) 



< he. 



□ 



As in [3], we need to divide the Lipschitz open set D into a good set and a bad set. We 
recall several geometric facts about the Lipschitz open set. 

Definition 4.5 Let e, r > 0. We say that G C 3D is (e, r)-good if for each point p G G, t/ie 
imzi inner normal v(p) exists and 

B(p, r) Pi (9.D C {x : |(x — p) ■ u(p)\ < e\x — p\}. 

If G is an (e, r)-good subset of then using this definition we can construct a good subset 
Q of the points near the boundary: 

g= \Jr r { P ,e), 



where T r (p, e) = {x : (x — p) ■ v(p) > y/1 — e 2 \x — p\} fl B(p, r). 

The next lemma is [3j Lemma 2.7] and it says the measure of the set of the bad points 
near the boundary is small. Note that even though [3j Lemma 2.7] is stated for a bounded 
Lipschitz domain, the proof remains true for a bounded Lipschitz open set. 

Lemma 4.6 (Lemma 2.7 in [3]) Suppose e G (0, 1/2), r > and that G is a measurable 
(e, r)-good subset of 3D. There exists so(dD, G) > such that for all s < s 

\{xeD: 5 D (x) < s} \ g\ < s [U d -\dD \ G) + e (3 + H d -\dD))] . 

The next lemma is about the existence of a good subset G C dD. Again the lemma 
remains true for a bounded Lipschitz open set D. 

Lemma 4.7 (Lemma 2.8 in [3]) For any e > there exists r > such that an (e, r)-good 
set G C dD exists and 

H d -\dD\G) <e. 
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The two lemmas above imply that 



\{x G D : 8 D {x) < s} \ Q\ < se (4 + H d -\dD)) . 



For any e G (0, 1/4), we fix the (e,r)-good set from Lemma [4.71 and construct Q from G. 
We choose r to be smaller than the minimal distances between (finitely many) components 
of D. For any x G Q, there exists p{x) G dD such that x G T r (p(x),e). Next we define inner 
and outer cones as follows 



I r (p(x)) = {y : (y-p(x)) -v(p{x)) > e\y - p(x)\} n B(p(x),r), 

U r (p(x)) = : (y -p(x)) -v{p{x)) < -s\y - p(x)\} (1 B(p(x),r). 
It follows from [31 (2.20)] that there exists a half-space H*(x) such that 

X G iT(x), <y ff . (a .)(x) = 5d(x), I r (p(x)) C F*(X) C U r (p(x)f . 



(4.8) 
(4.9) 

(4.10) 



Now we are ready to prove Theorem 11.21 
Proof of Theorem 11.21 Fix e G (0, 1/4), the (e,r)-good set from Lemma I4T71 and the Q 
constructed from G. From the definition of the trace we have 



-t d/a / r£(t,x,x)dx = t d/a / (p™{t,x,x) -p m (t,x,x))dx 
Jd Jd 

t d/a Z%(t) -t d/a [ p m (t,x,x)dx 
Jd 

t d/a Z£(t) -t d/a [ (p°{t,x,x) - (p°(t,x,x) -p m (t,x,x)))dx 
Jd 

t d/a Z 1 S{t)-C l \D\+t d/a [ (p°(t,x,x) -p m (t,x,x))dx. 

Jd 



Hence it follows from Lemma 13.21 that in order to prove Theorem 11.21 we must show that for 
given e G (0, 1/4) there exists a to > such that for any < t < t Q , 



f d/a I r ™( t ,x,x)dx - C 2 H d -\dD)t 1/a 
Jd 



< c(e)t 



l/a 



where c(e) — > as e — > 0. As in the proof of [3J Theorem 1.1.] we split the region of 
integration into three sets 

Vi = {x G D\g : 8 D (x) < s}, 

v 2 = {x g Dng : S D (x) < s}, 

V 3 = {xeD: 5 D (x) > s}, 
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where s must be smaller than the Sq given by Lemma 14.61 For small enough t we can take 

8 = t X/a /y/i. 

It is shown in [3j (3.2) and (3.4)] that 

t d/a [ r° D (t,x,x)dx < c(e)t l/a (4.11) 

where c(e) — > as e — > 0. Hence it follows from Lemma 12. II and (14. lip that 

t d/a [ r%(t,x,x)dx < c(e)e 2mt t 1/a . (4.12) 

Now we deal with the integral on V 2 . Let H*(x), I r (p(x)), U r {p(x)) be defined by (14 . 8 j) . 
gU) and (00]) . We have 

I r (p(x))cH*(x)cU r (p(x)) c . 
Since r is less than the minimal distances between components of D, we also have 

I r (p(x))cDcU r (p(x)) c . 
Since I r {p(x)) C U r (p(x)) c , By an argument similar to that used in Lemma [3.41 we have 

\rp(t,x,x) — rff*^(t, x, x)\ 

< e 2mt (rl ip[x)) (t,x,x) - r° Uripix)r (t,x,x)) . (4.13) 
Now it follows from [3, Proposition 3.1.] and (14. 13j) that 

^ ^ / I^ZJ 0*> "^J — ^ H* (x) (^'" 3>) | G^E 

</X>2 

< ce^'^Vv^^W 1 /" (r^ +1 M)dr. 
Finally we will show that the integral 

* d/a / r%., x) {t,x,x)dx 

gives the second term C2'H d ~ 1 (dD)t 1 / a plus an error term of order c(e)t l l a . Recall that 

r2. {x) {t,x,x) = f%*(t,5 H * ix) ) = f%(t,5 D (x)). 
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Hence we have 



JV 2 



t d/a [ f%(t,5 D (x))dx 



Jv 2 



t d/a [ f^(t,5 D (x))dx-t d ^ [ f™(t,6 D (x))d: 




By an argument similar to that used to get (14.121) we have that 




JT>iUT>3 

where c(e) — > as e — > 0. From the (approximate) scaling property of the relativistic stable 
process, we have 



Now apply Lemmas 14.21 and 14.41 to the function r — > /#*(!, r) and we get for small enough t 
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